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’ ABSTRACT 

The second-degree nonlinear aeroelastic equations of motion for 
a slender, flexible, nonuniform, Darrieus vertical-axis wind- turbine 
blade which is undergoing combined flatwise bending, edgewise bend- 
ing, torsion, and extension are developed using Hamilton's principle. 
The blade aerodynamic loading is obtained from strip theory based on 
a quasi-steady approximation of two-dimensional incompressible un- 
steady airfoil theory. The derivation of the equations has its 
basis in the geometric nonlinear theory of elasticity and the result- 
ing equations are consistent with the small deformation approximation 
in which the elongations and shears (and hence strains) are negligible 
compared to unity. These equations are suitable for studying vibra- 
tions, both static and dynamic aeroelastic instabilities, and dynamic 
response. Several possible methods of solution of the equations, 
which have periodic coefficients, are discussed. 


INTRODUCTION 

Renewed interest in the wind as an alternative source of energy 
has resulted in a number of studies of various wind-turbine concepts. 
Presently, receiving considerable attention (see for example Refs. 1 
and 2) is the vertical-axis wind turbine (VAWT), also known as the 
Darrieus rotor (Ref. 3). The VAWT rotor (Fig. 1) embodies long, 
slender, and flexible airfoil-shaped blades which are attached to a 
vertical rotating shaft at both ends. The curved shape of each blade’ 
is approximately that of a troposkien, i.e., the shape taken by a 
flexible cable of uniform density and cross section when it is spun 
at a constant angular velocity. 
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The efficient design, construction, and operation of large VAWT 
rotors require that the vibrating loads and stresses in the blades, 
as well as in the combined rotor-tower system be reduced to the 
lowest possible levels and that the system must be free from all 
types of instabilities. Thus, aeroelastic and structural dynamic 
considerations have a direct bearing on the manufacture, life, and 
operation of large VAWT systems. Although the basic dynamic phenomena 
associated with VAWT rotors are basically similar to those of helicop- 
ter rotors, proprotors, and horizontal-axis wind turbines, the struc- 
tural configurations of VAWT systems are sufficiently different to 
necessitate comprehensive and independent analytical and experimental 
investigations of their aeroelastic stability and dynamic response 
characteristics. 

Vertical-axis rotor systems can exhibit a variety of mechanical 
and aeroelastic instabilities such as resonance, ground resonance, 
whirl flutter, blade classical bending-torsion flutter, blade coupled 
bending-torsion-extension aeroelastic instabilities, blade stall 
flutter, and blade static divergence. Ground resonance and whirl 
flutter are associated with the entire VAWT system; the other insta- 
bilities are primarily associated with the individual blades of the 
system, A design requirement for VAWT systems is that each component 
as well as the entire system be free from all instabilities. 

Experimental results obtained with wind-tunnel models (Refs. 4 
and 5) have indicated that the blade aeroelastic instabilities involv- 
ing coupling between flatwise bending, edgewise bending, and torsion 
are possible under certain conditions. To explain these instabili- 
ties, an analytical investigation was conducted in Ref. 6 using an 
approximate modal analysis. More recently, bending vibration equa- 
tions of a rotating curved slender blade were derived and solved for 
special cases by using asymptotic methods in Ref. 7. An aeroelastic 
analysis of an existing 5m VAWT system with and without guying wires 
for the tower was conducted in Ref. 8. The analyses indicate the 
possibility of resonance, ground resonance, and aeroelastic-type 
instabilities. In Ref. 9, an analytical investigation of the aero- 
elastic stability of a different 5m VAWT system was performed using 
a finite element model. These results also show the possibility of 
several types of instabilities, depending on the system parameters. 

Analyses based on finite element models are well-suited to 
accommodate the structural complexities of actual VAWT systems which 
may have blades with struts and towers with guying wires. However, 
the fundamental understanding of the basic mechanisms of aeroelastic 
instabilities and dynamic response phenomena and parametric studies 
associated with VAWT systems are better served by a continuum model 
which leads to a set of differential equations. A continuum model 
for a single blade can be viewed as a building block from which a 
continuum model for an entire VAWT system can be constructed. 
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The purpose of this report is to develop a set of second-degree 
nonlinear aeroelastic equations of motion of a Darrieus wind-turbine 
blade involving flatwise bending, edgewise bending, torsion, and ex- 
tension. The nonlinear terms which will be retained in the present 
derivation are of the type which have been found to be important in 
aeroelastic stability of helicopter rotor blades. 

The derivation of the nonlinear equations of motion herein fol- 
lows the methodology of Refs. 10-12. The equations are derived using 
the geometric nonlinear theory of elasticity (Ref. 13) in which the 
elongations and shears (and hence strains) are negligible compared 
to unity. The generalized aerodynamic forces are obtained from strip 
theory based on a quasi-steady approximation of two-dimensional incom- 
pressible, unsteady airfoil theory. The equations of motion which are 
consistent with these approximations may be derived to any desired 
degree by retaining terms in the dependent variables to the appropri- 
ate degree throughout the development. The present development will 
be directed to the derivation of the second-degree nonlinear equations 
of motion in which one formally retains terms through second-degree in 
the dependent variables. Rigorous adherence to this retention scheme 
leads to an almost insurmountable amount of algebra. To circumvent 
this problem to some extent, an ordering scheme which is consistent 
with the assumption of a slender beam is imposed early in the develop- 
ment of the dynamic and elastic portions of the present equations. No 
ordering scheme is imposed in the development of the generalized aero- 
dynamic forces herein because any ordering scheme which is imposed 
would depend on the order assigned to the nondimensional free-stream 
velocity and induced velocity both of which vary significantly in 
practice. Thus, to accommodate such general operating conditions with 
the present equations, the aerodynamic forces are left in general 
second-degree form. The aerodynamic forces acting on a blade element 
are functions of the blade azimuth angle and hence the final equations 
will contain periodic terms. For completeness, the gravitational 
forces are also included in the present development. 

The equations developed herein are suitable for studying aero- 
elastic instabilities, aeroelastic response, and vibration character- 
istics of flexible, curved, and rotating blades. These equations 
form a building block from which a continuum model of an entire 
Darrieus-type VAWT system can be constructed. As these equations do 
not have closed form solutions, several possible approximate methods 
of solution are discussed. 


SYMBOLS 

a airfoil lift-curve-slope 

a-p bp c^ quantities defined in Eq. (A5) 


A 


cross-sectional area of blade 
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A T3 


A u’ ^v* \r 


B V’ B T’ B G’ b a 


B t (i = 1, 2, 3, 4) 


o 

C(k) 

D 


projected area of rotor in vertical plane 

generalized aerodynamic forces per unit length 
in ^x B 3 ’ ^yg3 # ^ z B 3 directions > respectively 

generalized aerodynamic moment per unit length 
about elastic axis 

number of blades 

boundary terms arising from strain energy., 
kinetic energy, work done by gravitational 
forces, and work done by aerodynamic forces, 
respectively 

sectional constants 

blade chord 

airfoil profile drag coefficient 

Theodorsen's circulation function 

airfoil drag per unit length 

chordwise offset of mass centroid from elastic 
axis (positive when in front of elastic axis) 

chordwise offset of area centroid of cross 
section from elastic axis (positive when in 
front of elastic axis) 

Young's modulus 


£X B2’ 

eY B2’ 6Z B2 

unit 

vectors 

along 

X B2 

» Y B2 ’ Z B2 

axes 

6x B3’ 

SY B3’ SZ B3 

unit 

vectors 

along 

X B3 

’ Y B3’ Z B3 

axes 

® X B6’ 

6Y B6’ SZ B6 

unit 

vectors 

along 

X B6 

’ Y B6» Z B6 

axes 



unit 

vectors 

along 

X I ’ 

Yj, Zj axes 

> 

£y 9 &<7 

y r z r 

unit 

vectors 

along 

X R 

Y^, axes 


aerodynamic force vector 
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F F F 
X B3* y B3* Z B3 


F F F 

X B 6 ’ y B 6 ’ Z B 6 


g 

G 


components of aerodynamic force vector F A in 
the directions of 'ey , ev , 'ey , 
respectively 


B3 


'B3 


J B3 


components of aerodynamic force vector F* in 


the directions of e 


X- 


B6 


ey 


B 6 


respectively 
gravitational acceleration vector 
shear modulus 


J B 6 


G u’ G V ’ G w 


gravitational forces per unit length in u, v, w 
directions, respectively 

generalized gravitational moment per unit length 
about elastic axis 


H 

h 


height of wind turbine 

vertical velocity of two-dimensional section 
normal to free-stream 


I , I , I 
U V w 


Ix 3 x 3’ Iy 3^3 


generalized inertia forces per unit length in 

, ey , ~e z directions, respectively 
B3 B3 B3 

area moments of inertia about Yg and 
axes, respectively 

generalized inertia moment per unit length 
about elastic axis 


'A 


torsional section constant 
reduced frequency 

polar radius of gyration of cross-sectional area 
about elastic axis 


k ± (i = 1 , 2 , ... , 6 ) 


k m 


^mj_ > k m? 


k , k , k 
X B3 y B3 Z B3 


notation used in writing the variation of the 
kinetic energy 

polar radius of gyration of cross-sectional mass 
about elastic axis (k^ = k^ + k^ 0 ) 


'•m2 > 


mass radii of gyration about Y 3 and X 3 axes, 
respectively 

components of curvature of elastic axis before 
deformation 
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k„ , k„ , k„ 
X B6 y B6 Z B6 


components of curvature of elastic axis after 
deformation 


li> tn-f 5 5 

(i = 1, 2, 3) 
L 


m 


M, M 


Z B6 


P 


P' 


r l 


R, Rl 


s, s x 

s i (1, 2,. ..,9) 



direction cosines, Eq. (Al) 

aerodynamic lift per unit length 

mass of the blade per unit length 

aerodynamic moment about the deformed elastic 
axis per unit length 

arbitrary point on the elastic axis before 
deformation; also origin of the blade-fixed 
axis system before deformation 

arbitrary point on the elastic axis after 
deformation 

position vectors of a point in the cross section 
of the blade before and after deformation, 
respectively 

position vectors of an arbitrary point on the 
elastic axis before and after deformation, 
respectively 

running coordinates along the elastic axis 
before and after deformation, respectively 

notation used in writing the variation of strain 
energy in a concise form 

length of blade along undeformed elastic axis 

generalized elastic forces 

time 


kinetic energy; blade tension; rotor thrust 


u, v, w 


U 

Up, Up, Up 


deformations of elastic axis in Xg 3 , Yg 3 , and 
Zg 3 directions, respectively 

resultant of Up and Up 

radial, tangential, and perpendicular components 
of the resultant velocity of a point on the 
elastic axis 
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induced velocity, positive in the negative Xj 
direction 


V 


strain energy 


V 


oo 


free-stream velocity 


V 


a 


wind velocity vector 


'y Y 7 ’ 

A B3 B3 B3 


V- 


X B 6 Y B 6 Z B 6 


relative velocity of a point on the elastic axis 
expressed in X B3 Y B3 Z B3 an< ^ X B 6 Y B 6 Z B 6 coordi- 
nate systems, respectively 


W 


sum of and Wq 


W 


A 


W 


G 


x 


o’ 


z 


o 


X 3 , y 3 , Z 3 


work done by aerodynamic forces 

work done by gravitational forces 

coordinates of a point on the undeformed elastic 
axis along X-j-- and Zj-axes, respectively 

coordinates in X 3 Y 3 Z 3 coordinate system, the X 3 
and y 3 axes are the minor and major principal 
axes of the cross section 


Wi 

X B1 Y B1 Z B1 

X B2 Y B2 Z B2 

X B3 Y B3 Z B3 

X B 6 Y B 6 Z B 6 


x r y r z r 




inertial axis system 

blade axis system, parallel to X^Y^Z^ coordinate 
system 

blade axis system obtained by rotating x bi y b1 z B 1 
system about the negative Y B1 ~axis by an angle 0 O 

blade axis system obtained by rotating x b 2 Y B2 Z B2 
system about the Z-g 2 ~ ax ^- s by an angle y 

blade axis system in the deformed configuration 

obtained by translating and rotating the x b 3 Y B3 Z B3 
system; the Zgg-axis is tangent to the deformed 
elastic axis 

blade axis system obtained by rotating the XjYjZj 
system about the Zj-axis by an angle \p(= fit) 

transformation matrix relating the angular orienta- 
tion of the deformed and undeformed blade-fixed 
coordinate systems 
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Green's strain tensor 


a 


airfoil section angle of attack 


a x , a^, a z quantities defined in Appendix A 

6, £, 9 Eulerian-type rotation angles 


Y 


section total pitch angle (built-in twist plus 
pitch angle due to control inputs) 


y x 3 z 3 , Y y 3 z 3 , Y z 3 z 3 

6 ( ) 


engineering strain components 
variation of ( ) 


69 virtual rotation about the Zp^-axis 

Z B6 B6 

e small parameter of the order of the bending 

slopes; airfoil section pitch angle with 
respect to free-stream velocity; also 
extensional component of Green's strain 
tensor along the elastic axis 


XoZ 


3 3 


£y 3 2 3’ 02-32 


3 2 3 


strain components 


angle between blade local tangent and vertical 
axis, illustrated in Fig. 1 


M 


nondimensional free-stream velocity, V^/^R 


nondimensional induced velocity, v^/ftR 


P 


mass density of the blade; also mass density 
of air 


°x 3 z 3 > *y 3 z 3 ’ 

P 

a>X B3 Y B3 Z B3 

0)y Y 7 

B6 B6 B6 


engineering stresses 

angle of twisting deformation about the elastic 
axis 

blade azimuth angle 

curvature vector of the undeformed elastic axis 
curvature vector of the deformed elastic axis 


fj 


rotational speed of rotor 
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^Ni 

(') 

! 

( ) 


circulatory aerodynamic term 

noncirculatory aerodynamic term 

time derivative ~ ( ) 
dt 

denotes differentiation with respect to s 


MATHEMATICAL MODEL AND ATTENDANT ORDERING SCHEME 

The mathematical model chosen in the present development is a 
continuum model. The presence of rotation introduces equilibrium 
centrifugal stresses which require the use of a geometric nonlinear 
theory of elasticity. There are several levels of approximation 
which may be considered in this theory (Refs. 10 and 13). The level 
of approximation used in the present development assumes that the 
elongations and shears (and hence strains) are negligible compared 
to unity. 

The wind- turbine blade considered in the present development 
consists of a slender, curved, nonuniform blade which can undergo 
combined flatwise bending, edgewise bending, torsion, and extension 
(axial deformation). The elastic axis, mass axis, and tension axis 
are taken to be noncoincident. The elastic axis is assumed to be 
coincident with the quarter-chord of the blade. The generalized 
aerodynamic forces are obtained from strip theory based on a quasi- 
steady approximation of two-dimensional, incompressible and unsteady 
airfoil theory. Gravitational forces are included. 

Blades presently being considered for VAWT applications have 
neither pretwist (built-in twist) nor control inputs for changing 
section pitch angle as do the blades for a horizontal-axis wind 
turbine. However, for completeness, a variable section pitch angle 
is included in deriving the second-degree expressions for the bending 
curvatures and twist rate and for the strains. 

An ordering scheme consistent with the assumption of a slender 
beam is introduced here to provide a systematic procedure for dis- 
carding higher-order terms while deriving the second-degree nonlinear 
aeroelastic equations. A mathematical ordering scheme was introduced 
in Ref. 11 for deriving the nonlinear equations for a slender heli- 
copter rotor blade. Considerations similar to those in Ref. 11 have 
been applied in the present report to establish an ordering scheme 
which is consistent with the slender curved blades of a VAWT rotor. 

In this scheme, a parameter e which is taken to be of the same order 
as the nondimens ional variables u/S, v/S, w/S, and cp is introduced. 
The order of the dependent variables and geometric quantities appearing 
in the equations of motion of this report are as follows: 
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u/S = 0(e) 
v/S = 0(e) 
w/S = 0(e) 
<fi = 0(e) 
a x = 0(e) 
a = 0(e) 

y 

a z = O(e^) 


x^/S = 0(e) 
y 3 /s = 0(e) 
z 3 /S = s/S = 0(1) 
0 o (s) = 0(1) 
x 0 (s) = 0(1) 
y 0 (s) = 0(1) 
z Q (s) = 0(1) 


( 1 ) 


By using the ordering assigned above, the order of the elastic 
and inertial terms which are retained in the second-degree nonlinear 
aeroelastic equations of motion of the curved blade considered herein 
are given in Table 1. The rationale for this scheme was discussed in 
Ref. 11. It should be noted that in the present development the ex- 
tensional deformation (w/S) is 0 (e) instead of 0 (e 2) as in Ref. 11 
because of the presence of the initial bending curvature. 


Freedom 

Extension 

Bending 

Torsion 


Table 1 - Ordering scheme 


Elastic forces 
0 (e 3 ) 

0(e 4 ) 

0 (e 5 ) 


Inertial forces 
0(e 3 ) 

0(e 2 ) 

0 (e 3 ) 


COORDINATE SYSTEMS AND MOTION VARIABLES 

Several orthogonal coordinate systems will be employed in the 
derivation of the equations of motion. Those which are common to 
both the dynamic and aerodynamic aspects are described in this 
section. 


1. Inertial system (I-system) Xj Y jZj — The Zj-axis of this 
system, as is shown in Fig. 1, coincides with the vertical axis of 
the shaft. The X^-axis is aligned with the free-stream velocity V^. 

2. Rotating system (R-system) XrYrZr — This system is obtained 

by rotating the I-system about the Z]--axis by an angle = fit, as 

is shown in Fig. 1. The shaft rotational speed is given by ft and 

is assumed constant. The coordinate transformation between the I- 
and R-sys terns is 
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j **li 

cos 

<P 

- sin i) 

0 

i ex R 

{ ^ = 

sin 

'P 

cos ip 

0 


Kj 

0 


0 

1 



3. Blade system 1 (Bl-system) Xg^Yg-^Zg^ — This local blade- 
fixed coordinate system, as shown in Fig. 1, is fixed to an arbitrary 
point, P, on the elastic axis of the blade. This frame translates 
along the blade elastic axis and it is parallel to R-system. 

% 

4. Blade system 2 (B2-system) Xb2 y B2Xb 2 — This system is ob- 
tained by rotating the Bl-system about the negative Yg^-axis by an 
angle 0 O , as shown in Fig. 1. The Xg£s Yg 9 , and Zg 2 axes are 

in the normal, binormal and tangential directions, respectively. The 
rotation angle 0 O can be obtained from the known geometry of the 
curved undeformed elastic axis in the XjZj plane by the parametric 
equation 

R = x o®X! + z o^Z x (3 > 


and is 


0 O = tan 1 (- x^/zq) 


(4) 


The coordinate transformation between the Bl- and B2-systems is 


e XBl 

/ — l _ 

\ e *Bl / ' 
I ® Z B1, 


COS 0 Q 
0 

sin 0 O 


0 

1 

0 


sin 0 £ 
0 

COS 0, 


; e X B 2 

\ ® Y B2 

%2 


(5) 


5. Blade principal axis system (B3-system) Xgg, Ygg, Zg3 — The 
Xg 3 and Yb 3 axes are taken to be aligned with the minor and major 
principal axes of the blade cross section, respectively. The prin- 
cipal axes are obtained by rotating the normal and binormal axes by 
an angle y as indicated in Fig. 2. The angle y is the total 
section pitch angle, which is a combination of built-in twist (pre- 
twist) and section pitch changes due to control inputs. The VAWT 
configurations presently considered in the literature do not have any 
section pitch angle, but it is included in developing the expressions 
for the curvatures and strains in order to indicate how one would 
include this effect in the analysis. The coordinate transformation 
between the B2- and B3-systems is 
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/ X 1 

r 

i 


“1 


A 

ex B2 

' cos y 

- sin y 

0 


Ex B3 

B2 

) 

i 

) = ; sin y 

cos Y 

0 

J 

\ 

^g 3 

1 

i e? 

' Z B2 1 

! 

0 

0 

1 


%3y 


The curvature vector of the undeformed blade is 


“ X B3 Y B3 Z B3 


- 0„ev + y ev 
° *B2 Z B3 


(6) 


(7) 


This vector is also called the Darboux vector in the literature. 
Substituting ey into Eq. (7) from Eq. (6), the curvature vector 
is given by 

wX Y Z = k v „e y + k„ ey + k„ e 7 (8 

B3 B3 B3 X B3 X B3 ^B3 Y B3 Z B3 Z B3 

where 


kx B3 


sin y 



0 Q cos y 


kz B3 


= Y 


! 


(9) 


6. Blade system 6 (B6-system) XggYggZgg — This system is shown 
in Fig. 3 and is obtained by translating and rotating the XggYggZgg 
system. The Zgg-axis is tangent to the deformed elastic axis. The 
blade cross section itself is assumed rigid. The deformations of the 
elastic axis are denoted by u, v, w in the B3-system. The angular 
orientation of the B6-system with respect to the B3-system is given 
by three Euler ian-type angles 8, and 0, which are, in turn, ex- 
pressed in terms of the elastic deformations (u, v, w) and their 
derivatives (u', v', w'), and the twisting deformation (<}>)• The 
final relation between the B3- and B6-systems is developed in 
Appendix A and is given by 


*S 

?X B6 

! _ ii _ 4 
1 2 2 




] 

- a x - <f>ot . 

j 

'N 

* X B3 

¥ y b6 / = 

- * - « x a y 

1 - 

2 

Uy 

T 

t 

\ 

(f>a x - CXy 

J 

j _ v 

\ %3 f 
i 1 

%6 ! 

a x 

a y 



i 1/2,2. 

1 - — (a + a ) 
2 x y • 

1 ! 

'. e 2 B3 1 


(10) 
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where 


a„ = u - vk„ + wk, T 
x Z B3 y B3 


c y - v + uk ZB3 - wk XB3 


a = w - uk + vk 
z y B3 B3 


(ID 


HAMILTON'S PRINCIPLE 


The governing equations of motion are derived using the extended 
Hamilton's principle (Ref. 14). 

t l 

(6T - 6V + 6W)dt = 0 (12) 

where 

6W = SW G + 6 W a (13) 

In Eq. (12), T Is the kinetic energy, V is the strain energy, and W 
is the work done by gravitational and aerodynamic forces. For subse- 
quent convenience, the variation of work is divided into two parts as 
indicated in Eq. (13): the first part, 6 Wq, is due to gravitational 

forces; and the second part, 6W A , is due to aerodynamic loading. In 
the following sections, explicit expressions for T, V, and W in 
terms of the dependent variables u, v, w, and 4> and their deriva- 
tives and the blade sectional properties will be developed. 



Strain Energy 

The expression for the strain energy of the blade in terms of 
engineering strains and stresses is 




r l 

JJ T z 3 z 3 Yz 3 z 3 
A 


+ 


a Y 

Z 3 X 3 ' Z 3 X 3 


+ cr 


z 3 y3 y z 3 y3) dx 3 dy 3 dz 3 


The 


coordinates 


s 


and 


Z3 are used interchangeably. 


( 14 ) 
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where, assuming that the components of the engineering strains are 
equal to the corresponding components of the Lagrangian strain and 
using Hooke's law, 


° 2 3 x 3 GYz 3 x 3 2 GEz 3 x 3 

0z 3 y 3 = GYz 3y3 = 2Gez 3y3 


(15) 


Taking the first variation of V as given in Eq. (14), and 
using Eq. (15) , yields 

I s r r 

6V = / E J J y Z3Z3 8 Vz 3Z3 d X3 d y3 dz 3 

o A 


+ 


J 


v 


(Y 


z 3 x 3^ Yz 3 x 3 + > dx 3'»'3< iz 3 


(16) 


The expressions for the required strain components are developed 
in Appendix B. For a slender curved blade with zero section pitch 
angle these expressions are given by Eqs. (B15), (B16) , and (B17). 
Substituting these expressions into Eq. (16), taking the indicated 
variations, and integrating over the cross section of the blade leads 
to 

n s 

5V = , (s]_<5u + S 2 < 5 u' + S 3 < 5 u" + s^ 6 v' 

(/ 0 


+ S 56 V" + Sgdw + Sy6w ' + Sg6(j) b Sgdtjl'jdzg (17) 


where, consistent with the ordering scheme discussed earlier, 
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So = T(u' - 9'w) 

* o 

S 3 - - e; [EEly 3 y 3 - 2EB 3 (v" + 4)0’)] 

+ Ei y 3 y 3 [u" - e>* - 0 ”w - ^ (V * 2 + (j) 2 ) + 4 >v"] 
+ T(f>e A - E4,T X3X3 (v" + 4>0q) + GJv' (4 1 - QqV' ) 


s 4 = Tv' + EIy 3 y 3 9iv' (u" - 9 qW - B^w') 

+ gj ( 4> ' - e^v')(u" - w'e^ - e^w - e^) - Gje^v’(u" - w’e^ - e^jw) 
s 5 = EI y y W - 0”W - 0>’) - Te A + 2EB 3 0 q (u" - 0>' - Bjjw) 


s 6 

s 7 

S 8 


- EI X 3 X 3 [~ V " + 4(u" - 0 o w - e i w ' - 0 o)] 

= Twe; 2 - T0’u’ - EI y3y3 0^(u" - 0^w - ey) + Ely^S'e 


= t 


+ 


+ 


-El 0 ’ (u" - 0"w - 9'w') 

^ 3^3 ° W ° ° 

Q ' 0 (v" + W[EIy 3 y 3 e - EB 3 (v" + 40’)] 

Ely 3 y 3 (v" + 0q4) (u" - 0Q W - 0>') + EB 3 (v" + 0i4) 2 0i 


E(u" - 9 q w ' “ 0oW)[Aee A - I X3 x 3 (v" + 49 q) 
0’B 3 ( u " - 9’w’ - e>)]- E0q {Ae A e - B 3 [0’(u" - 
Q o -f- - e’K'J + ^ Q’(4> 2 + 24V”) + Ix 3 x 3 [- V " 


0> - 0’w’) 

+ 4 (u" - 0 qW 


0 ’w 


') - 4>9o j + — 4' (4’' 

o j 2 


20oV')} - 2B 3 E0^(u" - 0 qW - 0>') 2 


(cont ' d) 
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- sA 2 e {-b 3 


u" - Qqw' - 6qW + <f>v" + j ©o (4> 2 + v' 2 ) 


Srs = 


- V 0 o v " + e o 2 ^j + E6 ^ v " + 6 >>j I x 3 x 3 e - V v " + K> 

(4>* - e;v')[E(Iy 3 y 3 + Ix 3 x 3 ) £ - E(B 3 + B 4 )(v" + <J>6q)~] 

+ GJjj) ’ + v'(u" - w' 0 q - e^w - e^)’j 


T * E iAe - I 


y 3 y 3 


6' 2 <j> 2 

9o(u" - 9> - 0QW') + -+ e^v" 


+ — ~ 0^(0^ + 2v" ) + AeJ-v" + <Ku" - 0> - 0^w' - 0^)J 
2 

i 

+ ~r t'(*' - 20'v’) 

2 (18) 


The expression for the extensional strain e on the elastic 
axis for the case in which the section pitch angle is zero is given 
by 


e = w’ + Q^u + -r (u' 2 + 0^ 2 w 2 - 2 GqU'w + v' 2 ) (19) 

Assuming that the cross section is symmetric about the Yg 3 ~axis, the 
sectional properties appearing in Eq. (18) are defined as follows: 


A ■/ / dx 3 d y 3 



x 3 dx 3 d Y 3 


0 


Ae A = ,/ ‘f y3 dx 3 d y3 


x 3 y 3 dx 3 dy 3 = 0 


I 


X 3 X 3 



y 3 dx 3 dy 3 



(cont 5 d) 
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y 3 y 3 


Ak" 


o 0 


x^dxgdy^ 


U 0 


Pi o 2' 

/ H + yi idx 3 dy 3 



y|) dx 3 d y 3 


>"3 ( x 3 + y 3) dx 3 dy ; 


y^x^dx^dy^ = 0 


r) n 

x y^dx dy = 0 
i ./ . I 3 3 3 3 


u , , y 3 x 3 dx 3 d!, 3 ■ 0 


b 3 - 



x^y^dxgdy^ 


n 


B 4° J L / y 3 dX 3 dy 3 


( 20 ) 


Integrating Eq. (17) by parts, the result can be put in the form 

n s 


6V = 


(S u 6u + S v t 5 v + S w 6w + S ( j ) 6<fi)dz3 + By 

Jo 


(21) 


where the generalized elastic forces S u , S y , S w , and S. are 


S u S 1 s 2 + 8 3 


S v = ~ s 4 + S 5 


J w *6 “ s 7 
S <J) = s 8 ~ s 9 

and the boundary term is given by 


( 22 ) 


= 1 ( s 2 - S3)5 u + S36U' + (S4 - scpdv + s^dv' + sydw 4 - sgdcf> j 


0 

( 23 ) 
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Kinetic Energy 

The position vector of an arbitrary mass point of the blade is 


r 


1 


R + AR + x~e v + y~e v 

J X B6 J Y B6 


( 24 ) 


Using Eqs. (5), (10), and (A16) in Eq. (24), the expression for r^ 
with respect to the B3-system is 


where 


r i X 3®X. 


B3 


+ Ye + Z e 

0 133 1 ^ B 3 


(25) 


/ 2 

x 3 = x D cos 0 O + z 0 sin 0 O + u + X3 (l ~ 2 ~ ~ y ) ~ y 3^ + a x a y) 
<j> 2 a y\ 

Y 3 “ v + y 3 1 - f + *3$ 


= - x Q sin 0 O + z D cos 0 O + w - x 3 (a x + <j) 0 Sy) + y 3 (<j>a x - ay) 

(26) 


The angular velocity of the B3-system can be written as 


u = d,e 7 - fie 7 = £!e 7 (2 7) 

L I ^R ^B1 

Substituting for e from Eqs. (5) and (6), yields 

2 B1 

a) = Ci sin 0_ cos ye Y - sin 0. sin ye v + cos 6 n e 7 j (28) 
o a b3 o x b3 0 ^B 3 / 

The section pitch angle y is set to zero in the subsequent develop- 
ment. For this special case Eq. (28) simplifies to 


go = ii sin 6 n e Y + fi cos 0 o e 7 (29) 

° X B3 0 Z, B3 

_ The expression for the kinetic energy of the blade in terms of 
r^ and u is given by 
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T 



dr l 

— dx 3 dy 3 dz3 


(30) 


where 



r-j_ + a) x r]_ 


(31) 


The variation of T, integrated between t Q and t^, is given by 


{ 



1 


6T 




dx 3 dy 3 dz 3 dt 


(32) 


Substituting Eqs. (25) and (29) into Eq. (31) and the result into 
Eq. (32), integrating by parts over time where necessary, and then 
integrating over the cross section, the variation of T can be put 
into the form 


ns 


6T = / (k^6u + k 2 <Su' 4* k 3 6v + k^dv' + k 3 dw + kgd<j>)dz 3 

^0 


(33) 


where, consistent with the ordering scheme introduced earlier, 

.. ... 2 2 
k^ = - m(u - e<j>) + 2m ftv cos 0 Q + mfi cos ® 0 ( X D cos 0 O + u - e $) 

2 

- mfi sin 0 O cos 0 Q (-x sin 0 O + w - ev') 

2 2 2 2 

k^ = - mv’ex Q fi cos 0 Q + mk^ft sin 0 O cos 0 Q - mfre<J>x 0 sin 0 Q 

,, m # • * o 

k 3 = - mv + 2mfi sin 0 o (w-ev') - 2mfi cos 0 Q (u ~ e<J>) + mSF (v + e) 

2 ... 2 
k^ = - me(u' - 0 qW)Q x q cos 0 q + me(w + 2vft sin 0 Q ) + mfi ex Q sin 0 Q 

70 9 

- ms ew sin* 1 0 Q + mil eu sin 0 O cos 0 Q 


(cont ' d) 
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= m Ik^^G^Cu' - 0 qW ) ft 2 cos^ 0 Q + v'0q [ -cu + 2evft cos 0 Q 

2 2 2 i 

+ eft cos ® Q Cx o cos 0 Q + u) - eft sin 0 Q cos 0 o (-x o sin 0 Q + w)j 

- w - ecf>(u' ~ 0 qW) ~ 2e4>(u' - 0 qW) + ev' - ej>(u' - S^w) 

2 9 r 

- 2ft sin 9 o (v - e<j>4> ev'v') + ft sin z 0 Q j -x Q sin 0 q + w 

+ e<f>(u' - 0 qw) - evQ - ft 2 sin 0, o cos 0 O | x 0 cos 0 O + u 

f 2 .. 2 . 

- ev'(u' - 0 qw) - etj) j + 0 q k (u' - 6 qw) - 2ftk^ (p sin 0 

ill ^ O 

- ft 2 sin 2 0 q k^Cu' - e^w) - ft 2 sin 0 Q cos 0 Q k^ 

- 0q | - ew - 2ftev sin 6 q + ft 2 sin 2 0 O e(-x 0 sin 0 Q + w) 

2 

- ft sin 0 O cos 0 Q e(jj D cos 0 O + n 

2" 2 2 9 

Kg = meu - - 2mftev cos 0 Q - mft ex Q cos e Q - Ml eu cos^ 0 Q 

+ mft 2 ew cos 9 Q sin 0 Q + m k^ 2 - k^jtjjft 2 cos 2 0 O 

m-i (^2 ~ ^2.^ V S ^ n & ° C ° S 0 ° ~ m (^2 ~ ^ 

-- n , .eft 2 x 0 (u 1 - 6 qW) sin 0 Q -- 2mk 2 ^ft(u' - 0^w) sin 0 Q 

-.2 . •• , 

- me;i ‘(j)v + 2meftdm cos 0 O + me<f>v - 2meftw<f> sin 0 Q 

- mew(u' - SqW) - 2meftv sin 0 o (u' - 0 qw) - meft 2 u sin 0 Q cos 0 o (u' 

- 9'w) + meft 2 w sin 2 6 0 (u' - O^w) 




( 34 ) 
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The sectional properties appearing in Eq. (34) are defined as follows 


'/-< 


m 


pdx 3 dy 3 me = J py 3 dx 3 dy 3 


wk l 1 = JJ P x 3dx 3 dy 3 mk 2 


m2 U ( 




py 3 dx 3 dy 3 


•4 ■ > 4 , + 4 2 


'P 

I px 3 y 3 dx 3 dy 3 = 0 


pdx 3 dy 3 = 0 


V 


(35) 


Integrating Eq. (33) by parts, the resulting expression can be put in 
the form 


6T - ./ (I u 6u + I v <$v + I w 6w + I 6 5<J>)dz 3 + B T 


(36) 


where 


lu " k i - k 2 


I v = k 3 - ki 


X w = k 5 


= k 6 


and the boundary term is given by 


(37) 


Brj. = (k2<$u + k^Sv) 


(38) 


Virtual Work of Gravity Forces 

The virtual work due to gravity can be expressed in the form 

1 S 11 n 

pf • 6t r 1 dx 3 dy 3 dz 3 (39) 


6W g - 


0 A 
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where g is the gravitational acceleration vector and is given by 


g = -ge Z][ (40) 

The vector "g can be expressed with respect to the B, 3-system by 
substituting for e Z j from Eqs. (5) and (6). For the special case 
of zero section pitch angle considered in the present development, 
the vector "g in the B3-system simplifies to 

S " " 8 ( ¥ X B3 sin 9 P + ®Z B3 GOS 0 ; < 41 ) 


Taking the variation of r^ which is given in Eqs. (25) and (26), 
substituting the resultant expression together with Eq. (41) into 
Eq. (39) , integrating over the cross section, and integrating the 
result by parts, Eq. (39) yields 


6W & — j (G u 6u + G v 6v + G w 6w + G ! j ) 6i)dz3 + Bq 


(42) 


where 


G u = - mg sin 6 Q + (- mgev' sin 0 O + mge# cos 0 Q ) 

I 

G v = |- mge(u' - 0qw) sin 06 - mge cos 9 0 ] 


G w = - mgev ' 0 q sin 0 O - mg cos 0 O + mgeO^ cos 0 Q 


C, = mge sin 0 - mge cos 0_ (u* - 0 'w) 


(43) 


and the boundary term Bq is given by 


B 


q sin 0 Q - mge# cos 0 Q ) 6u + jjnge(u' - 9 q w ) sin 0 Q 

S 


+ mge cos 0 Q <5v 


0 


(44) 
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Virtual Work of Aerodynamic Forces 


The virtual work of the aerodynamic forces can be written as 





+ M 60 
Z B6 



dz 3 


( 45 ) 


where 6Y-^ is the virtual displacement of the position vector of an 
arbitrary point on the elastic axis, is the aerodynamic force 

vector, 60 z gg is the virtual_rotation about the Zg^-axis. Usually 
the aerodynamic force vector F^ is calculated in the B-6 blade 
axis system. Since the position vector Y^ given by Eq. (25) is 
expressed with respect to the B3-system, the force vector F^ is 
transformed to the B3-system using the following relation 


iW 


X B3 Y B3 Z B3 


= M {Fa> 


X B6 Y B6 Z B6 


(46) 


The aerodynamic force in the Zgg .direction is Fzgg and is a pro- 
file drag force. Following usual practice, this force component is 
assumed to be unimportant and is taken to be zero. Substituting 
Eq. (A39) into Eq. (46) and discarding terms which will lead to terms 
higher than second-degree in the final equations, one obtains 









+ F 


y B6 


, « - F_ (u’ - elw) - F v' 

Z B3 X B6 0 y B6 


(47) 


Taking the variation of the position vector r^ (Eq. (25)) on the 
elastic axis yields 


6r-i = 6ue Y + 6ve y + 6we 7 (48) 

1 X B3 B3 X B3 

The virtual rotation 60 Z gg is obtained from the expression for 
k Z gg given in Eq. (A40) by replacing Qq by 60 o , cf> ' by 6<j>, and 
(u' - 0 q w ) by 6(u' - 6 qW) , and making use of the fact 60 o equals 
zero, and is 
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60 = 6(f) + v ! (6u' - 6'6w) 

Z B6 0 


( 49 ) 


Substituting Eqs. (47-49) into Eq. (45), the virtual work expression 
reduces to the form 


6W A = 


(A u 6 u + Ay.dv + A^w + A ( j ) 6 <J))dz 3 + B A 


0 


where the generalized aerodynamic forces are 


(50) 


A = F -F d> - (M v') 
u *B6 YB6 Z B6 


A = F + F d> 
v y B6 X B6 


^ = - F x b6 (u ' - - F Y B6 V ' " 


A , = M 
* Z B6 

and the boundary term is 


B a = M v'6u 
A Z B6 


(51) 


(52) 


There remains the task of expressing 




and 



in 


terms of the dependent variables u, v, w, <j> and the geometric angle 
0 Q . These expressions will be generated from two-dimensional, incom- 
pressible, quasi-steady, strip theory in which only the velocity com- 
ponents perpendicular to the span-wise axis (zgg-axis) of the deformed 
blade are assumed to influence the aerodynamic loading. Account will 
be taken of the pulsating free-stream velocity V(t) associated with 
a rotating blade by employing Greenberg's extension of Theodorsen's 
unsteady theory (Ref. 15) for determining the aerodynamic lift and 
pitching moment acting on the blade. The resulting expressions are 
specialized to the case of quasi-steady flow by setting Theodorsen's 
circulation function to unity. Classical blade element momentum 
theory is used to calculate the steady flow induced by the rotor. 


In the present application of Greenberg's theory, the airfoil is 
taken to be pivoted in pitch about the aerodynamic center at the 
quarter chord and to be executing harmonic motions in pitch (e(t)) 
and plunge (h(t)) while immersed in a pulsating airstream V(t), as 
shown in Fig. 4. The lift and moment acting on the elemental section 
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of the blade may be expressed in terms of the circulatory and non- 
circulatory components as 


L = L c + 

M = M c + Mnc 


(53) 


Since the blade elastic axis is assumed to be coincident with the aero- 
dynamic center at the quarter chord, the individual components of 
Eq. (53) follow from Ref. (15) and can be written as 


NC 


M, 


‘NC 


M, 


~ pa ^h + Ve + Ve + ej 
pacV ^h + Ve + | ej 

- j (ff + '^ + fs) 

- i pac (ff 2 Ve 


(54) 


In the course of arriving at the circulatory terms in Eq. (54), the 
quasi-steady approximation has been introduced by setting the reduced 
frequency to zero, in consequence of which Theodorsen's circulation 
function C(k) assumes the value of unity. 


The lifts and moments given in Eq. (54) must now be expressed in 
terms of Ur, Up, and Up, the radial, tangential, and perpendicular 
velocity components relative to a point on the elastic axis of the 
airfoil, Fig. 5. Now, the expression in parentheses for L^c in 
Eq. (54) is the downward acceleration relative to the air of the mid- 
chord point of the airfoil, and the expression in parentheses for Lq 
is the downward velocity relative to the air of the three-quarter- 
chord point of the airfoil. Since Up is the relative velocity com- 
ponent perpendicular to the quarter chord, the sectional lifts can 
also be written as 


L 


NC 


A ■ * c •• 

2 pa 4 i-Up + j e 


2 pacU (-Up + 2 £ 


(55) 
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where V(t)» appearing outside the parentheses of the expression for 
Lq in Eq. (54), has been approximated by the resultant of only the 
tangential and perpendicular velocity components and is given by 


V - U = (56) 

As indicated in Fig. 6, the noncirculatory lift acts normal to the 
section chordline and the circulatory lift acts normal to the result- 
ant velocity U. The profile drag force acts parallel to U and is 
given by 

1 ^dn 9 

D = j pac U 2 (57) 

where is the (constant) profile drag coefficient. 

u o 

The components of the aerodynamic force in the directions of the 
X B6’ Y B6’ and Z B6 axes are given by 


F x = + Lq cos a - D sin a 


YB6 


1^ sin a - D cos a 


F Zb6 = “ D r (neglected) 


where 


(58) 


sin a = Up/U 
cos a = U T /U 

Substituting Eqs. (55), (57), and (59) into Eq. (58), leads to 


(59) 


F *B6 ' i PaC i “ U f D T + I V - f °P + (t 

F y B6 - I oac i U P -f V -T°1 0 ! 




UU T 


(60) 


The noncirculatory and circulatory moments in Eq. (54) can be written 
in terms of Up, Up, U, and z and assume the form 
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( 61 ) 


from which the total pitching moment is given by 

2 

M = M = - pac (y) (lie - U_ + — ce) (62) 

zg6 2 V 4 / \ P 8 I 

The necessary expressions for Up, Up and e are developed in 
Appendix C, and are given by Eqs. (C14) and (C19) . Substituting 
Eqs . (C14) and (C19) into Eqs. (60) and (62) and the resulting ex- 
pressions into Eq. (51), one obtains the necessary expressions for 
generalized aerodynamic forces. 


AEROELASTIC EQUATIONS 


Expressions for (ST, <5V, and <SW have been obtained in the pre- 
vious sections. Substituting these expressions and their associated 
boundary terms into Eq. (12), there results an expression of the form 



£( ) 6 u + ( )<5v + ( ) 6 w + ( ) 6 <jT]dz 3 + ® 7 dt = 0 


(63) 


For arbitrary admissible variations, 6 u, 5v, 6 w, and <5cj), the four 
expressions in parentheses must vanish individually as must the 
assembly of boundary terms denoted by B. The first condition will 
yield the four governing nonlinear partial differential equations 
for u, v, w, and <f>, and the second condition will give the associ- 
ated boundary conditions at the ends of the blade. The governing 
equations of motion and boundary conditions are summarized below. 


u equation: 


m(u - e<£>) - 2 mfi cos 0 q v - mfi^ cos 0 o x o - mfi 2 cos^ 


2 -22 
+ mfi sin 0 cos 0 (w - ev') +1 mk n sin 0 O 

O O j 


2 . o 

- men v x n cos 0 - meQ^x sin 0 d> 


+ mg sin 


0 o (u - ecj>) 
cos 0 

o 

0 Q + (mg ev' sin 0 Q 


(cont'd) 



2EB 3 (v" + <f>e£) j + EI y y 


eA 


u" - 0Aw’ - 0"w - (V - 2 + <j> 2 ) + <pv" I 


o 2 

M 


+ T<j»e A - E4>I X3X3 (v" + + GJv ' (<j> 1 - 0£v') J. = ^ 


v equation: 

••2 

mv - mfl (v + e) - 2mfi sin 0 Q (& - ev') + 2mfi cos 0 o (u - e4>) 

+ [mew + 2mevfi sin 0 Q + mn^ex Q sin 0 O - mfi 2 ew sin 2 0 Q 

T 

+ mfi^eu sin 0 Q cos 0 O - mefi 2 (u' - 0 Aw)x q cos 0^] + |jnge cos 0 O 

T 

+ mge(u’ - 0;w) sin 0 Q ] - [Tv* + ' (u" - 0^w - 0^w') ' 

+ GJ(<fc' - 0^v' ) (u" - w'0q - 0qW - 0^) - GJ0Av'(u" - w'©A - 0^w)] 

+ | EI y y 4>(u" - e'Aw - e^w') - Te A + 2EB 3 0A(u" - - e^w) 

- ei x 3 x 3 '--v" + ♦(«" - 0> - ey - ©AflJ = \ 

w equation: 

mw + me<j)(u' - 0 qw) + 2me^>(u' - 0 qw) - mev' + meij)(u’ - 0 qw) 

2 £ 

+ 2mQ sin 0 o (v - ecjxp - ev'v') - mfi sin 0 Q |-x D sin 0 + w 
+ e<Ku' - ©A w ) - ev^J + mft 2 sin 0 Q cos 0 Q [x Q cos 0 q + u - ev'(u' 

- Q' q w) - ecj>] - m^^eACu' - ©Aw)fi 2 cos" jj-eu+2ev^ cos 0 q 


(cont ' d) 
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+ eft 2 x Q cos 0 O + efi 2 u cos 2 d Q - eQ 2 w sin 6 o cos Qq] 


+ m 0 Q 


2 A ' V 2 A,' - ft 't 


-k 2 (u ' - 0>) + 2Bl4 1 * sin 0 O + n? sin z 0^ k^Cu' - 0^w) 


+ Q 2 kj^ sin 0 Q cos 0^_ 


+ mcj)0Q j^-ew - 2ftve sin 0 Q - fi 2 ex 0 sin 0 Q 


+ n 2 sin 2 0 ew - £2 2 sin 0_ cos 6 eu ; + Tw0' - T6'u' 
o o o — u u 


- ei VoVo 0 q(u" - 0 > - fl>') + Eiy^^o’-ei- 


-y 3 y 3 


- eAw') 


T ” EI y 3 y 3 0 ° (u " “ e ° w 


+ mg cos 0 Q + mge sin SqV'Oq - mge©^ cos 0 0 = ^ 


4> equation: 

•« • 2 2 2 
mk 2 d) - meii + 2mev cos 0_ + mefi x n cos 0 + mef^ u cos 0 


m 


me2 2 w sin 0 Q cos 6 0 - m \*l 2 ~ ^lj ^ ^ °° + \ k m 2 


->2 ! 


/ 2 


- kj^Jv' sin 0 Q cos 0 Q + m fkJ 2 - kJ^B 2 # + mefi^x 0 sin 0 o (u' 


0'w) 

o 


+ 2mk^L ft sin 0 Q (u' - O^w) + meft 4>v - 2meftcf>u cos 0 O - me<j>v 
+ 2meftwi}> sin 0 Q + mew(u' - 0 gW) + 2meftv(u' - S^w) sin 0 O 


+ meft 2 u(u' - 0 qW> sin 0 O cos 0 Q - meft 2 w(u' - 0 qW) sin 2 0 Q 


e ;(v- + 0 ^) 


0> - 0’w') 


EI y 3 y 3 E " eb 3 (v " + * e o>j + EI y 3 y 3 (v " + e i* ) 

+ EB 3 (v" + Kr,V) 2 ^' Q + E(u" - 6 ^w - e^»')[Ae A E 


u 


Ix 3 x 3 Cv" + *9') - 6'B 3 (u" - e>’ - 0»] - E0' / Ae A e 


(cont ’d) 
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0Q 2 (j> 2 

9 o(u" - 9 > - 9 >') -“V - - 6 >’ 


B i 9 n 9 

+ -y 2 (0 2 + 2<j>v") 


+ I xoX9 l_- v " + 0 (“" - e'ow - eiw') - <f>ei] + ^ - 20^') 


X 3 X 3 


,2 j 


2BoE0'(u" - 0"w - e'w') z - E0' z { -Bo (u" - 0'w' - 0"w + 0v" 
Jo v o o' o|3| v o o T 


+ \ 0^(0 2 + v' 2 ) 


- B 2 6q(v" + 0q 4>) ) + E0^(v" + e^4>) 


I £ 
X 3 X 3 


- B 4 (v" + ©o4>) | - { (0* - 9 o v ' ) j E (l y ^y 3 + 1x3X3) E " E(B3 + B4)(V " 


+ 0 9 o) 


i I 


+ GJ <j>' + v' (u" - w'0q - 0q'w - Q ' q ) ! \ - mge sin 0 


+ mge(u' - 6 qW) cos 0 = 
where 

e = w' + 0^u + j (u' - 0 qW) 2 + | v' 2 


(64) 


X = E (ae - I e i(u" - 0 qW - e>’) + 0' 2 £ + 0’0V" 


*3*3 


+ I x 3 x 2 ( e o <t> + 2v ") + Ae AL _v " + 0( u " " e o w 


0 n W ’ 

o 0 


- e’l 

oJ 


+ Ak 2 ^ (0* - 20'v')j 


(65) 


It should be noted that when Eq. (65) is substituted into Eq. (64) 
some third-degree terms in u, v, w, and 0 and their derivatives 
result. Since only second-degree equations are of interest herein, 
the third-degree terms should be discarded. Similarly, when Eqs. (C14) 
and (C19) are substituted into Eqs. (60) and (62) and the resulting 
equations into Eq. (51) some terms higher than second-degree in u, 
v, w, and 0 and their derivatives result. These terms should also 
be discarded in the final expressions for A u , A v , A w , and A^. 

The assembled collection of boundary terms denoted by B is 
given by 


B = B-p - By + Bq + B a 


( 66 ) 
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and the requirement of the vanishing of the individual variational 
components leads to the relations 


(s 1 - S3) - k 2 - (mgev’ sin 0 O - mge<f> cos 0 O ) - M z v 


6u 


s^du' 


= 0 


(34 - S5) - k 4 

,S 
0 


s^dv' i = 0 


mge(u' - 0'w) sin 0 n + mge cos 6 f 


dv 


= 0 


St6w = 0 

'O 


sgd<j> 


= 0 


(67) 


METHODS OF SOLUTION 

The governing equations of motion are coupled, nonlinear, partial 
differential equations with periodic coefficients In the dependent 
variables u, v, w, and <f>. These equations have no closed-form solu- 
tion and must be solved using approximate methods. Usual practice in 
solving these equations is to, first, eliminate the spatial dependence 
This results in a set of coupled nonlinear ordinary differential equa- 
tions with periodic coefficients. Various techniques can then be 
employed to solve these equations to determine either aeroelastic 
stability or response. Some of these techniques are briefly summari- 
zed below. The reader interested in detailed considerations should 
refer to the references cited. 

The spatial dependence is usually eliminated by employing a 
modal approach with either assumed or calculated mode shapes (Ref. 14) 
An alternative procedure for eliminating the spatial dependence Is by 
use of an integrating matrix approach (Refs. 16-18). 

The nonlinear ordinary differential equations with periodic coef- 
ficients also have no closed-form solution and must be solved by 
approximate methods. A common practice has been to numerically inte- 
grate these equations in time to determine time histories of u, v, w, 
and 4 > from which aeroelastic stability and response of the system 
can be determined. The assessment of stability can be facilitated if 
fast Fourier transforms are performed on the time histories (Ref. 19). 
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Blade moments and shears can also be calculated using the time his- 
tories. Another method for solving the nonlinear ordinary differential 
equations has been the classical perturbation method for determining 
stability whereby the nonlinear equations are perturbed about a steady- 
state equilibrium position. This leads to two sets of equations: a 

set of nonlinear algebraic equations for the steady-state quantities 
and a set of nonlinear ordinary differential equations with periodic 
coefficients in the perturbation quantities. The nonlinear algebraic 
equations are solved by standard iterative techniques. Usual practice 
is to linearize the perturbation equations by discarding all the per- 
turbation terms of second-degree or higher in the perturbation vari- 
ables. These linearized equations are then solved using Floquet- 
Liapunov theory (Refs. 20-22) from which aeroelastic stability and 
response can be determined. Another method sometimes employed for 
solving the linear perturbation equations is an approximate solution 
based on time-averaged coefficients in conjunction with a so-called 
multi-blade coordinate transformation (Ref. 22). In this approach, 
the linear perturbation equations are first transformed into a non- 
rotating coordinate system by means of a multi-blade coordinate trans- 
formation where some of the periodicity in the coefficients is trans- 
formed into constant terms. A constant coefficient approximation is 
then made by time averaging the remaining periodic coefficients in 
the differential equations. Standard eigensolution techniques can 
then be used to determine aeroelastic stability and response. 


CONCLUDING REMARKS 

The second-degree nonlinear aeroelastic equations of motion for 
a slender, flexible, curved, and nonuniform Darrieus vertical-axis 
wind-turbine blade undergoing combined flatwise bending, edgewise 
bending, torsion, and extension have been derived using the extended 
Hamilton's principle. The blade aerodynamic loading is obtained from 
strip theory based on a quasi-steady approximation of two-dimensional, 
incompressible, unsteady airfoil theory. The derivation of the equa- 
tions has its basis in the geometric nonlinear theory of elasticity 
and the resulting equations are consistent with the small deformation 
approximation in which the elongations and shears (and hence strains) 
are negligible compared to unity. A mathematical ordering scheme 
which is consistent with the assumption of a slender beam was adopted 
for the purpose of systematically discarding higher-order terms in the 
elastic and dynamic forces in the final equations of motion. The 
expressions for the generalized aerodynamic forces were left in gen- 
eral second-degree form from which one can obtain the aerodynamic 
loading to the order appropriate to any case of interest. The final 
equations, which have periodic coefficients, are suitable for studying 
vibrations, both linear and nonlinear aeroelastic stability and 
response. As these equations do not have closed form solutions, 
several approximate methods of solution have been discussed. 
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APPENDIX A 


DERIVATION OF EXPRESSIONS FOR COORDINATE 


TRANSFORMATION MATRIX AND CURVATURES 

The development of the second-degree nonlinear aeroelastic equa- 
tions of a rotating curved blade requires second-degree nonlinear 
expressions for the rotational transformation matrix between the B3- 
and B6-systems (Fig. 3) and for the components of curvature. Since 
these expressions are independent of rotational speed, ft, only a non- 
rotating blade is considered. 

The elastic deformations translate and rotate the B3-system to 
the B6-system. Let the rotational transformation matrix between these 
two systems be [t] such that 


fey ^ 

x Bb 


' ex B3' 


*1 

m l 

n l 

r— *\ 

SX B3 

* Y Bb 

> - M < 

B3 

>- 

*2 

m 2 

n 2 

\ 6 y B3 
1 



^- ez B 3> 


_f 3 

m 3 

n 3 

i 

t 

1 67 

Z B3 ^ 


Let the expression for the curvature vector of the deformed elastic 
axis be 


W X 


B6 Y B6 Z B6 


= k 


' x B6 SX B6 + ^Bb^Bb + ^Bb^Bb 


(A2) 


The next step is to find the expressions for the components of 
the curvature vector in terms of the direction cosines mi, 

np ... n 3 , and the components of initial curvature. From Eq. (Al) , 
one can write 


e Z_ = *3 e X_ + m 3 e Y„ „ + n 3 e Z 


J Bb 


B3 


'B3 


B3 


Differentiating Eq. (A3) with respect to s leads to 


J Bb 


= a 2 e X +bie v f c,e 


B3 


l e Y 


B3 


1 Z 


B3 


(A3) 


(A4) 


where 
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a-t = - m-jk, + n^k 

'B3 


1 ~ X, q — UioJtV„ '* liqK.,, 

1 3 3 z DO 3 y B 3 


bi = + & 3 k z __ - n-^k. 


B3 


3 K x- 


B3 


C 1 " n 3 + m 3 k Xg 3 " ^3^y 


'B3 


(A5) 


Invoking the small strain assumption, one can write 


der. 


J B6 


ds- 


= (1 + 2s) 


de Z 


ds 


B6 

" e-. 


J B6 


(A6) 


The identity 


leads to 


-17 OJ-y V 7 X ^-7 

Z B6 X B6*B6 Z B6 z B6 


(A7) 


SZ B6 ky B6 6X B6 kx B6 £Y B6 


(A8) 


Substituting Eq. (Al) into Eq. (A4) , one can also write e' in the 
form 


e' = (a,H, + b-im-i + c,n-,)e 
Z B6 1 1 X B6 

+ (a-^ + ^l m 2 + c l n 2^®Y 

B6 


+ (a^£ 3 + b^m^ + c^n-^e^ 


B6 


(A9) 


From Eqs . (A8) and (A9) , the expressions for k Xg6 and k yi!C are 


TB6 


kx B6 = " (a l £ 2 + b l m 2 + Cin2) 


ky B6 = 3l£l + bimi + Cini 


(A10) 

(All) 


The identities 
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C-y- — U)y y rj ^ 6y 

a B6 *B6 B6 4 B6 a B6 


K e 7 + K e Y 
y B6 Z B6 Z B6 B6 


Z B6 ” X B6 " Y B6 


(A12) 


and the substitution of Eq. (Al) into Eq. (A12) lead to 


? 

k = l^i’ev + m-ilv + nie 7 | • i&oev + m 0 ev + n-Vey | 

Z B6 V 1 X B3 1 X B3 1 Z B3 / V z X B3 1 X B3 Z ^B3j 


(A13) 

Expanding Eq. (A13) , using Eq. (8), and making use of the fact that 
each element of the orthogonal matrix [t] is equal to its cofactor, 
Eq. (A13) simplifies to 


J B6 


£ 3 k x + m 3 k v + n 3 k 7 
B3 y B3 B3 


+ ^ 2^1 + m 2 m l + n 2 n l (A14) 


Thus far, the expressions for the components of curvature have been 
developed in terms of the direction cosines and the components of the 
initial curvature of the undeformed elastic axis. The next task is 
to express the direction cosines in terms of u, v, w, and <j>. To 
this end, the direction cosines are first expressed in terms of the 


Eulerian-type angles, 8, £, ■ 

and 

0 which are defined as follows: 

1. A positive rotation 
X B4 Y B4 Z B4’ 

8 

about 

the 

Y B3 -axis 

resulting in 

2. A positive rotation 
ing in X B5 Y B5 Z B5 . 

C 

about 

the 

negative 

Xg^-axis result- 

3. A positive rotation 
X B6 Y B6 Z B6‘ 

0 

about 

the 

Z B 5~axis 

resulting in 


The explicit form of the transformation matrix [t] in terms of the 
Eulerian-type angles is 
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cos 0 cos 3 

sin 0 cos £ 

-cos 0 sin 3 

-sin 0 sin £ sin 3 


-sin 0 sin £ cos 3 

-sin 0 cos 3 

cos 0 cos £ 

sin 0 sin 3 

-cos 0 sin £ sin 3 


-cos 0 sin £ cos 3 

cos £ sin 3 

sin £ 

cos £ cos 3 


(A15) 

The quantities 3, £, and 9 are to be expressed in terms of 
the variables u, v, w, <j>, and the components_of the initial curva- 
ture of the elastic axis. To this end, let R and be the 

position vectors of the points P and P T in Fig. 3, and AR be 
the displacement of P. Then, 


R-i = R + AR = R t ue v + ve v + we 7 (A16) 

1 X B3 *B3 ^B3 

Differentiating the expression for AR in Eq. (A16) with respect to 
s leads to v 


^ + «y«Y B3 + “z e Z B3 


(A17) 


where 


a„ 

= u' 

- k„ v 

+ 

k, T w 

X 

Z B3 


?B3 

a, y 

= v' 

+ k„ u 


k„_ w 

y 


Z B3 


X B3 

a z 

= w ' 


+ 

k Y „ v 
X B3 

(A16) 

with 

respect 

to s 


(A18) 


Eq. (A17) into the resulting expression, the expression for SR-^/Ss 
is 


"1 SR 


So - „ + ct„ey + ot„ e 7 

Ss Ss x X B3 y * B3 z Z B3 


(A19) 


Now, from calculus 


SR = - 
Ss e Z B3 


(A20) 
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Substitution of Eq. (A20) into Eq. (A19) gives 


3R, 


- u, c > a e v + (1 + a )e„ 

3s x X B3 y^ B3 z Z B3 


= a e„ + 


The relation between the extensional component of the Green 
strain tensor on the elastic axis and SR^/Ss is given by 


(A21) 

’s 


E = 


1 mi 3Ri 


3s 


2 \ 3s 


- 1 


(A22) 


Substituting Eq. (A21) into Eq. (A22) , the expression for e reduces 
to 


e = a z + \ ( a z + a x + a y) (A23) 

The relation between ds and ds^ can be written as 


1/2 

ds x = (1 + 2e) ds 


(A24) 


Differentiating Eq. (A16) with respect to s-^, and substituting 
Eq. (A24) into the resultant expression, one can write 


dR ± 

3sj_ Sz B6 


- 1/2 


(1 + 2e) 


(1 + a_ 


) e 7 + 

Z B3 


a„e v + 
x X B3 


x y e Y 


B3j 


(A25) 


Invoking the assumption that the elongations and shears (and hence 
strains) are negligible compared to unity leads to 


J B6 


“ a x e X.o + a y e Y B3 + (1 + a z> e Z 


B3 


B3 


From Eqs. (Al) , (A15) , and (A26) , one can write 

& 3 = cos 4 sin B = a x 
m^ = sin 4 = 

n 3 = cos ? cos 3 ** 1 '+ a 


(A26) 


(A27) 
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The orthogonality condition between £ 3 , m 3 , and n 3 

2 2 2 

^ + 112 = 1 (A 28) 

is satisfied before invoking the small strain assumption. This con- 
dition must also be satisfied under the small strain assumption. 
However, Eq. (A26) which assumes small strains leads to an interesting 
result. From Eqs. (A23) and (A27) , one can write for small strains 

n 3 = l + a z =l+e--|- + a y + a 2 ) - 1 - ^a 2 + a y + a z) 

(A29) 

An alternative expression for 03 in terms of a x and oty follows 
from Eqs. (A27) and (A28) and is 


(•* + «?) 


n - 


(A30) 


Thus, Eqs. (A29) and (A30) show that there are two slightly different 
expressions for 03 . These two expressions must be equal. Therefore, 


one should impose an additional assumption that the quantity at is 


negligible compared to 
assumption is invoked, 
ment. Accordingly, 


and and/or unity when the small strain 

This assumption is made in the present develop- 


= 1 - 


a 2 + 
x 


a 2 ) 

y! 


(A33) 


From Eqs. (A27) and (A31), the expressions for the trigonometric func- 
tions involving B and 5 are 


sin 8 = x — - a„ cos 3 = 1 - ~7r 

1 - Oty/ 2 X 2 

2 

0t\7 

sin 5 = a^. cos 5 = 1-—- (A32) 

The retention of the terms a x /2 and a^/2 in Eqs. (A31) and (A32) 
is consistent with the fact that some of the rotations must be regarded 
as substantially exceeding the strain components for a slender rotating 
beam. This implies that the right hand side of Eq. (A23) represents a 
small difference of large terms. This is discussed in Ref. 13, 
page 203. The implications of discarding these terms while deriving 
the nonlinear aeroelastic equations of a helicopter rotor blade were 
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discussed in Ref. 11. The third Eulerian angle 0 is associated with 
torsion of the blade and, hence, is given by 


0 = <j> 


(A33) 


The expressions for the transformation matrix [tJ and for the 
components of curvature are given in terms of the direction cosines 
in Eqs . (A10) , (All), and (A14) , and those for direction cosines in 
terms of the Eulerian angles are given in Eqs. (Al) and (A15) . The 
Eulerian angles are expressed in terms of the quantities a x , ay, a z , 
and <f> in Eqs. (A32) and (A33) , and the expressions for a x , a y , and 
a z are given in terms of u, v, w, kxg^, ky B g, and k Z gg in Eq. (A18) . 

Combining Eqs. (Al) , (A15) , (A32) , and (A33) , the second-degree expres- 
sion for the transformation matrix is 


[T] = 


(f, 2 a x 

2 2 

<f> 

-a x - <J>a y 


2 „2 


- 4> - a x a y 

! _ <*» _ ^ 
2 2 

* a x - a y 

a x 

“y 

i 

CM X\ 

s 1 

r— 1 



2 2 J 


(A34) 


Combining Eqs. (A5) , (A10) , (All), (A14) , and (A34) , the second-degree 
expressions for the components of curvature are 


k = k “ a x k 7 " a v “ k x ~ k x + ^ a x “ ^ a v- k z , + 

X B6 X B3 B3 y X B3 2 X B3 2 x y y 


B3 


B3 
(A35) 


(a 2 + ip 2 ) 

k = k - k a + - k — 2— L + <f>a’ + k Ja x 

y B6 y B3 B3 7 y B3 2 y Z B3 


Z B6 


k a x + k a v + k - k 
X B3 x y B3 y Z B3 Z B3 


lCy C6y C6 tt 

Y X B3 X B3 x y 

(A36) 

a y) 

+ <P' + a y a x 

(A37) 
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For the case of zero section pitch angle, Eq. (9) leads to 

k v = - 9A 

y B 3 ° 


k„ = k_ =0 
^B3 Z B3 


and Eqs. (A34) and (A35)-(A37) simplify to 


<j> 2 (u' - 0 qw) 2 


- (u* - 0>) 


w 


- (u' - 9 qW ) V 1 


(u' - 0lw) 


1-4-^r 


V 


1 - 


(u’ - e 


k x B6 = - V" + <J>(V - 0»' - <()0o 


’ 0 O , ,2 ^ ,2, 


yB6 


' + (u' - e'w) + -o- (v' + <f> ) + <j>v' 


Z B6 


0qv' + <p' + v' (u' - 8 qw) 


(A38) 


- 4>v * 


’w) 2 + v' 2 | 
(A39) 


(A40) 
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APPENDIX B 


DERIVATION OF STRAIN-DISPLACEMENT RELATIONS 

This Appendix will develop the second-degree nonlinear expres- 
sions for the strains. To this end, let r Q and r-| be the position 
vectors before and after deformation of an arbitrary mass point in the 
cross section of the blade. These vectors can be written as 


r„ = R + 


x 3 e X B3 + y3 e Y 


B3 


r ! " R ! + x 3 e X B6 + y3 e Y fi6 


(Bl) 


where the effect of warping of the blade is neglected. 
■ The differential of the vector, r Q , is given by 


dR 


de 


X 


B3 


dey [ 
Y B3 


dr o =[ d¥ + x 3 “d r~ + y 3 ~^j ds + e x B3 dx 3 + %, d y 3 (B3) 


B3 


The derivatives of the unit vectors can be written as 


de- 


X- 


B3 


ds 


a)x B3 Y B3 Z B3 Sx B3 


(B4) 


de-v 


ds “ X B3 Y B3 Z B3 * SY B3 


(B5) 


Substituting Eq. (8) into Eqs. (B4) and (B5), and the resulting equa- 
tions into Eq. (B3) , yield 

d7 0 = ( dx 3 " y 3 k z B3 ds > e X B3 + ^3 + x 3 k z B3 ds > e Y B3 


+ (1 - Xnk + y,k )dse„ (B6) 

j y B3 o x B3 ^ B3 


The same procedure leads to the following expression for dr 1 
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dr. 


dx~ - 


y3 k 2B6 ds l)% 6 + ( dy 3 


x 3 k z 


B6 


ds 


O' 


C B6 


+ (1 - x 3 k + y 3 k )ds,e (B7) 

' J y B6 J X B6/ 1 Z b6 


The Lagrangian 


strain tensor 



is defined as follows: 


dx 3 " 

dr l • dr x - dr 0 • dr Q = 2 [_dx 3 dy 3 dz 3 J [e^] ' dy 3 > 

ld Z3 J 


(B8) 


Combining Eqs. (A24), (A34) , (B6) , (B7) and (B8) and collecting terms, 
the expressions for the three strain components of interest become 


Z 3 Z 3 


- x„ I k 

- k 

\ + yofk 

- k 

3 ' y B6 
J2 

y B 3 . 

) 3 v X B6 
2 

X B 3 

^(k J - 

k 2 ^ 

<N 
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1 

+ 

- k 2 

2 \ y 36 

y B3 / 

2 V X B6 

X B 3 


A* A) 


Z B6 Z B 3 


x„y „ (k k - k k \ 
3 3 V y B6 X B3 y B3 ) 


(B9) 




-z 3 x 3 2 ’ Z B6 Z B3 


(BIO) 


£z 3 y 3 = ^ ' kz B 6 " S 3 


(Bll) 


Substituting Eqs. (A2 3 ), (A35), (A 3 6), and (A37) into Eqs. (B9) , (BIO) 
and (Bll) , the second-degree expressions for the strain components 
reduce to 


2 r 

** 3*3 *°« + i(“5 + 4) + f - 2k y B 3 k *B 3 “y + " k y B3 * 2 


+ 2k da' + 2k^ k„ „<j>a v - 2k„ k, r <f> + 2k Y k„ a„4> 

yb 3 y tb 3 zb 3 x x b 3 yB 3 x b 3 z b 3 y 


I v 2 I k 2 *2 _ 


2 ,2 


2k x B3 K + i yj ! ^ 2 - 2k x R , k z^ a x - 2k x RR “; - W 


2 J I yB 3 


X B 3 z B 3 


c B3 y X B 3 
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+ 2k k <f> - 2(j>a x k k - 2k k 4»a + 2k <fra^ - 2k <f>a ' 
X B3 y B3 y B3 Z B3 B3 Z B3 y B3 y B3 y 


xj- k a y + a’ - 

3 L Z B3 2 


yB3 ^a 2 + <J> 2 ) + <f>a' + k <t>a - <f>k. 
' y > y Z B3 x 


*B3 


ky, a a 
-x b3 x y 


+ 


y 3 


- k a - a* - — k^ - k + <f>a' 

Z B3 x y 2 ^B3 2 ^3 x 


< pa k_ + d>k„ 

y Z B3 y B3 


- x 


3 y 3 


- k k a v + k a’ - k k 

B3 B3 y B3 x B3 y B3 2 


+ 


2 

k v <j>a' + k k 7 <t>a - k„ <j> - a k k + V k a <p 
X B3 7 B3 B3 x ^3 y B3 B3 ^33 B3 


2 2 
+ k„ ^cba^ _ k„ „k v „„ §_ k„ da' - k„ k„, <J>a, + k <f> 

^B3 y X B3 y X B3 y B3 2 + yB3 x y B3 253 y y B3 


+ <f,k - k k <l>‘ 

y B3 B3 y B3 J 


IA 1 A) 


4> ' 3 + 2k„ k„ a 
X B3 Z B3 x 


+ 2k a*' + 21c, k a + 2k a <j»' + 2k <j) ' 
x b 3 x ^r^ z rr y ”rr y Yr^ 


S3 z B 3 y y B3 yY yB3^ 


(B12) 


y 3! 


- 23 X 3 


k„ a„ + k^ a. 

y-R'i 


2 |_ x B3 X y g3 y z B3 


2 4. 

a Y + a v J , ,1 

— . — ZJLL + <j) ' + a a ' 

2 y x 


(B13) 


X' 


Ez 3 y 3 2 i kx RR ax + ^R^y k; 


B3 


B3 y Z B3 2 


K + °y l +t > 


A ' + a a ' 


(B14) 


It should be pointed out that in arriving at the expressions given in 
Eqs. (B12) to (B14) several terms have been discarded based either on 
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considerations related to the small strain assumption or on consider- 
ations related to the approximations which can be made because of the 
assumed slenderness of the blade. 

For the case of zero section pitch angle, the expressions for 
the initial curvature components are given by Eq. (A38). Substi- 
tuting Eqs. (A18) and (A38) into Eqs. (B12)-(B14), the expressions 
for the strains are 


ZoZ 


= w' + e^u + \ (u' 2 + 0^ 2 w 2 - 20 q U ' w + v' 2 ) 


+ -T - 2e o (u " " 6 o w - 0 o w '> - Q o V - 20'<j)V r 


2 

Y3 

+ T ( 0 o V + 26'0 v") - x 3 


u" - 0^w - e>’+^ (V* 2 + (j) 2 ) + <j>v" I 


+ y 3 |- v" + cj>(u" - 0 qW - 0qW') - 4>0^ 


- x 


3 y 3 


- 20'cf,(u" - 0"w - 0; w ') + 8 q 2 (| 


2 , 2 


0 ; V " 


^ X 3 + ^3 ) 

+ 2 ■ ' (<S>' - 20^’v') 


(B15) 


'Z3X3 


5 [c’ + v’Cu* 




0 oW) - 0 O V '1 


(B16) 


’3^3 


= ~T + v '( u " - w’0Q - 0 qW> - e .v0 


(B17) 
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APPENDIX C 


DERIVATION OF THE VELOCITY COMPONENTS OF A BLADE ELEMENT 

The resultant velocities of a point on the elastic axis of the 
blade in the deformed and the undeformed coordinate systems are related 
according to 


? X Y Z * M'x Y Z 
B6 B6 B6 A B3 B3 B3 


where, from Fig. 5, 


X_Y “ U P e X_ , U T®Y_ + U R e Z 


B6 B6 B6 


B6 B6 


B6 


(Cl) 


(C2) 


and [t] is given in Appendix A. The total relative velocity (aero- 
dynamic + dynamic) of a point on the elastic axis is given by 


'Y Y 7 

A B3 B3 B3 


— dr-, 
V a - 1 


dt 


-l X B3 Y B3 Z B3 


(C3) 


Neglecting wind shear and gusts in the wind, the aerodynamic velocity, 
V a , consists of two parts: (1) the free-stream velocity V^; and 

(2) the induced velocity v^. As shown in Fig. 1, both and -Vj[ 

are parallel to the Xj-axis. Thus, 


^^XjY-j-Z-j- 


= (V 


v i )e "x. 


Defining two nondimensional parameters 


(C4) 


v = v /a R (C 5 ) 

CO 

and 

= v-^/ftR (C6) 

Eq. (C4) reduces to 


^ Y a^v y 7 fiR(p 

Vri 


p i^ e X I 


(C7) 
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Substituting for e^- in the above equation from Eqs. (2), (5), and 
(6), Eq. (C 7 ) reduces to 


(V ) = (y - V^ftR 

a X Y Z 1 

B3 B3 B3 


(cos 0 Q cos ip cos y - sin y sin if))e. 


X. 


B3 


- (cos 0 cos if) sin y + sin if) cos y)e v - sin 0 cos ipe? 

B3 ° Z B3 

For zero section pitch angle, Eq. (C8) simplifies to 


(C8) 


(V«). 


X B3 Y B3 Z B3 


= (p - y.)ftRjcos 0 cos tfie - sin if) e- 


B3 


B3 


- sin 0 cos ihe„ 
o Z 


B3 


(C9) 


For zero section pitch angle, the position vector of a point on 
the quarter-chord point from Eq. (25) is 


^K 3 m y 3 =0 ' <X ° C ° S ®° + Z ° Sln e » + “ )e X B3 + Ve I B3 


+ (- x 0 sin 0 O + Z 0 cos 0 Q + w)< 


J B3 


(CIO) 


The angular velocity of the B3-system from Eq. (29) is 


to = ft sin 0 e v + ft cos 0 e, 


^B3 


oZ T 


(Cll) 


The dynamic velocity of a point on the quarter-chord point, from Eqs. 
(CIO) and (Cll), is 


d^l 

dt 


- (u - vft cos 0 o )e„ + j_v + (x Q sin 0 O - z D cos 0 O - w)ft sin 0 


B3 


+ (x 0 cos 0 0 + Zq sin e o + u)ft cos 0 O ] e v + (w 4- vft sin 0 o )e 7 

- *B3 


(C12) 
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X B3 Y B3 Z B3 


■g Eqs. 

(C9) 

and (C12) 

into Eq 

. (C3) yields 

_ 1 
'L 

(y - 

%) 

OR cos 8 

0 

COS if) - 

(u - vft cos i 

+ 

!_- 

(y - 

y^)0R sin 

if) - v - 

(x„ cos 0_ + 
v o 0 

+ 

u)n 

cos 

8 o - < X C 

sin 0 - 

0 

z cos 0 - 1 

0 0 

- 

o 

- y. 

l> sin 9 0 

COS if) + 

(w + vO) sin 0 ( 


,)]■ 


c3 


J B3 


I B3 

(C13) 


Substituting Eqs. (A39) and (C13) into Eq. (Cl) and using Eq. (C2) , 
the second-degree expressions for Up and Up in terms of u, v, 
w, and $ and their time derivatives, are 


Up = - (y - y^)OR cos 0 q cos if) 


2 (u - 0 o w) 

1 - ±- - 

2 2 


2 


+ (u - vO cos 0 Q ) + <()(y - y^) OR sin if) + v<j> + x o Q<j) - w(j>0 sin 
+ dun cos 0 Q - L(u' - QqW) + <j>v'3 (y - yp)nR sin 0 Q cos if) 

- (u 1 - 0 q w) (w + vO sin 0 O ) 

U T = v'(u' - 0’w)] (y - yp)fiR cos 0 Q cos if) - <f>(u ~ vO cos 0 Q ) 


2 2 \ 2 1 2 \ 

+ (l - - ~-J (y - y^JiR sin if) + v + ;1 - - ~-j x q 0 


- wO sin 0„ + uO cos 0 + (y - y.)OR sin 9 cos if) (- v' 

O O 1 O *- 

+ <j)(u' - 0 q w )3 _ v' (w + vO) sin 0 O (C14) 

The quantity e appearing in Eqs. (60) and (62) is the angular 

velocity of the blade section about the local negative Zgg-axis and, 

consistent with the present notation, can be written as -e, . The 

Z B6 



48 


expression for £zb 6 can re § ar( ied as composed of two parts: the 

first part arising from the angular velocity of the shaft in space; 
the second part arising from the angular velocity associated with the 
elastic deformations. The first part is due to ft and can be obtained 
for zero section pitch angle from the relation 


f e x 

X B6 


f il sin O.'') 
o 

e yB6 

!> - M \ 

0 

v E ZB6> 


ft cos 0 
V o J 


(C15) 


Substituting Eq. (A39) into Eq. (C15) , the first part of e z is 

*B6 


(% 6 L- <u ’ - 


0'w)ft sin 0_ + 
o o 


1 _ (u’ - e;w)' 


,t2 


ft cos 0 


(C16) 

The second part is due to elastic deformation and is obtained by re- 
placing <j>' by (u' - 9qw) ' by -J— (u'-e^w), and 0^ by t) 0 (in 
the first term only) in the expression for k Z gg given in Eq. (A40) . 
Since 0 O is zero, the second part simplifies to 


\ Sz B6, 


deformation 


+ v’ (u' - 0lw) 


(C17) 


Combining Eqs. (C16) and (C17) , the total sectional pitching velocity 
of the airfoil is 




- 0qW)O sin 0 Q + 



(u* 


- e^w ) 2 
2 


v 


,2 


2 | 


ft cos 0 C 


+ (f> + .V ' (u' - 0 qW) 


(C19) 


Thus far, the expressions for Up, U T and e have been developed 
in terms of the quantities u, v, w, <j>, and their derivatives, the free- 
stream velocity the rotational speed ft, the geometric properties 

of the blade, and the induced velocity v-j_. Now, an explicit expression 
for the induced velocity will be developed. 
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The induced velocity Vp is determined from momentum theory in 
which the induced velocity at the rotor is one half its value in the 
wake. With this assumption, the thrust of the VAWT is 

T = 2pA p (V ro - v i )v i (C20) 


where Ap is the projected area of the rotor in the vertical plane. 
Another expression for the thrust will be developed in terms of the 
elemental forces acting on a blade section. 

Using Eqs . (2), (5), (47), and (60), the elemental force acting 
on a blade section in the Xp directions is 


% = ( F *B6 " Se*) cos * cos e ° ■ (% 6 * + F yJ sin * 


+ 



(u’ - e'w) + F v v’ 
° y B6 


sin 0 O cos ip 


(C21) 


The expression for the average thrust can be written as 



(C22) 


Equating Eqs. (C20) and (C22) , one obtains an integral equation for 
Vp which can be solved by an iterative procedure for a given V ro . 
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Figure L - Vertical-axis wind turbine and coordinate systems. 
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Figure 2. - Coordinate systems of blade cross section. 
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